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Abstract 



We show that the dynamic approach to Levy statistics is characterized by 

aging and multifractahty, induced by an ultra-slow transition to anomalous 

scaling. We argue that these aspects make it a protoptype of complex systems. 
PACS numbers: 05.40.fb; 05.45.Df; 05.45.Tp; 89.75.Da 
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Typeset using REVT^ 



The dynamic approach to Levy statistics is still a poorly understood problem in spite of 



from Continuous Time Random Walk (CTRW) |Q to the Nakajima-Zwanzig projection 
method This letter establishes that the dynamic approach implies aging, multiscaling 
and multifractality, thereby explaining, among other things, why this issue escaped so far 
a satisfactory understanding. The aim of this letter is to contribute to the comprehension 
of this delicate issue. Before addressing the problem from a more technical perspective that 
will be by necessity less accessible to a general audience, we want to illustrate the problem, 
and the solution of it afforded by this letter as well, with intuitive and qualitative arguments. 
To make our illustration as clear as possible, we adopt the same tutorial approach as that 
used by the authors working in the field of Cryptography |^, and we introduce two new 
archetypal individuals. Bob and Jerry. We hope that Bob and Jerry might have, with aging 
induced multifractality, the same fortune as Alice and Bob with Cryptography. These two 
individuals aim at realizing a process of diffusion of Levy type, with an intensity that depends 
on secrete numbers known only to them. Since the width of the probability distribution does 
not depend only on time, but also on the intensity of the jumps made by the random walkers 
at any unit time, it is hard, in principle, to establish the time at which the diffusion process 
began, if the intensities of the jumps are not known. However, as we shall see, this is 
possible with Bob's experiment, while it is impossible with Jerry's experiment. In fact. 
Bob and Jerry realize Levy diffusion in two different ways, and, as we shall see. Bob, who 
adopts the dynamic approach, allows us to predict the exact time at which he started the 
experiment, while it is not possible to guess the right starting time in the case of Jerry's 
experiment. Bob generates a random sequence of pairs {rj, Si}, with i = 0,1, ...oo. The first 
number of each pair, Xj, is randomly drawn from the distribution 



The index /x has to fit the condition jj, > 2, which ensures that the mean waiting time 



several attempts made at deriving it from intermittency with several techniques ranging 




(1) 
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is finite. Thus the number T controls the intensity of (r). The second symbol, Sj, is a sign, 
+ or — , and it is obtained by tossing a fair coin. Let us imagine that this distribution is used 
to generate a diffusion process according to the following prescription. Bob, who has at his 
diposal a virtual infinite number of walkers, creates, for any of them, a sequence {rj, Si} and 
makes her travel with velocity SiW for the time r^. Then Bob selects a number r belonging 
to the interval (0, tq]. The space travelled by the walker at time t > tq — r is given by 

x{t) = W[{to - r)so + Tisi + • • • 

(3) 

+Tn-iSn-i + (ttv - r')sN], 
where r' = (tq + • • ■ + r^v) — r — t and denotes the number of drawings made by Bob for 

his walker within the time interval [0, t]. At times t < tq — r {N = 0) the space is given by 

x{t) = WiTo - r - r')so, (4) 

with r' = Tq — r — t, yielding a;(0) = 0. Bob keeps secret both the value of W and the value 
of T. He can try to make his diffusion process look older by increasing either both or only 
one of these secrete numbers. Furthermore, to erase any possible form of aging, he selects 
the number r randomly. In fact, this has the effect of ensuring the stationary condition used 
by the authors of Refs |jl|J^ to realize Levy diffusion, under the form of Levy walk [|l| , which 
seems to be more realistic than the flight prescription 0,0. 

To appreciate the properties of the Levy walk, realized by Bob's experiment, it is con- 
venient to contrast it with Jerry's experiment. Also Jerry has at his disposal a virtually 
infinite number of random walkers, whose position at t = is x = 0, and he too, for any 
of his random walkers, selects an infinite sequence {xi}. The numbers Xi are drawn from a 
symmetric distribution, the positive numbers having the same probability as the negative 
numbers. For this reason Jerry does need the coin tossing to select Sj. Furthermore, Jerry 
makes his walker jump at any time step, by a jump of intensity \xi\, in the positive or nega- 
tive direction according to the sign of Xj. To be more precise, let us say that the distribution 
used by Jerry is Il{x), defined through its Fourier transform. 



= exp(-6|A;|^-i), (5) 

and only Jerry knows the secrete value of b. The distribution of numbers used by Jerry is 
the well known Levy distribution |^ : a stable distribution yielding a diffusion process, with 
the probability distribution pL{x,t), whose Fourier transform is 

PLik,t) = expi-b\kr~H). (6) 

The time t is the number of drawings, but it is so large as to be virtually indistiguishable from 
a continuous number. It is clear that the observation of Jerry's diffusion process does not 
allow any observer to establish when he began his experiments. We assume that the observer, 
which might be a third archetypal individual, does not know the time at which Jerry began 
his diffusion. By means of the experimental observation he/she can only establish bt, and 
since b is not known to him/her, he/she cannot determine the value of t. In other words, a 
broad distribution can be the consequence of Jerry starting his diffusion process at a very 
early time, but it can also be the consequence of a late beginning with much more intense 
jumps. 

It is not so with Bob's diffusion process. Let us see why. Let us consider, as in the case of 
Jerry's experiment, a time t very large. In the case of interest here, /i > 2, the mean waiting 
time, Eq. (j^), is finite. Thus, the number of random drawings and coin tossings is very 
well approximated by = ^/(t). Using the Generalized Central Limit Theorem (GCLT) 
we predict that Jerry's experiment yields the same statistics as Bob's experiment. Levy 
statistics. However, this important theorem does not afford any clear indication about the 
time necessary to realize this statistics. The predictions of the GCLT theorems are realized 
T0| by the following expression for p{x,t): 



p{x,t) = K{t)pT{x,t)e{Wt - \x\) + ^6{\x\ - Wt)Ip{t), (7) 
where 9 denotes the Heaviside step function. We also note that 

lim/,(t) = $^(t)= f-^)' \ (8) 



t-oo ' «^ ^ \T + t 



PT{x,t) is a distribution that for t oo becomes identical to the anti- Fourier transform of 
Eq. and K{t) is a time-dependent factor ensuring the normahzation of the distribution 

We observe that at any time t, no matter how arbitrarily large, it is possible to find a 
significant number of Bolb's walkers, with tq — r > t, namely, walkers for which Bob has not 
yet drawn the second pair of stochastic numbers. This probability is expressed by |[Tl|] 



It is interesting to notice that, due to the fact that Bob decides the motion direction by 
tossing a coin, the function ^^(t) is the correlation function of the variable C,{t), namely, the 
fluctuating velocity created by Bob's experiment. The number of walkers contributing to the 
propagation front is slightly larger. However, it is straightforward to prove with arguments 
similar to those used by Zumofen and Klafter that in the asymptotic time limit Ip{t) 
becomes identical to thereby accounting for the former of the two equalities of Eq. 

(|). The latter is easily accounted for by plugging iP{t) of Eq. (|I]) into Eq. (^. 

On the basis of these arguments we reach the conclusion that in the asymptotic time 
limit Eq. (|^) becomes identical to (see for earlier derivation) 



This equation, although valid only in the asymptotic time limit, is very convenient for the 
theoretical arguments of this letter. First of all, it allows us to determine the age of the 
diffusion experiment created by Bob, even if Bob adopts the stationary condition and 
keeps secret the values of W and T. To do so, we measure the distance of one ballistic 



two ballistic peaks. All these three quantities can be expressed in terms of the unknown 
quantities, t, W and T. The distance between the two peaks is: 2Wt] the diffusion coefficient 



is given by: b = W{TWY'hin[{TT{fi - 2)/2]r(3 - jj,)) ^ and the peak intensity by Eq. 
(I). The age of Bob's experiment can be revealed by means of experimental observation. 




(9) 



p{x, t) = pl{x, t)e{Wt -\x\) + -6{\x\ - Wt)<!>^{t). 



(10) 



peak from the other, the diffusion coefficient h of Pl(x, t) of Eq. (^ and the intensity of the 



5 



thanks to the breakdown of homogeneity (multi-scahng), caused by the ultra-slow relaxation 
of At the intuitive level of this first part of the letter, we can say that aging causes 

inhomogeneity. 

Let us now move to a more technical level of description. We share the vision of Khinchin 
]r3| about the close connection between ordinary statistical mechanics and the ordinary 
Central Limit Theorem. There is, on the other hand, a close connection between statistical 
equilibrium and scaling of a diffusion process. The latter property is expressed by 

pix,t) = ^F{^). (11) 

In fact, this property means that the probability density at different times can be expressed 
in terms of the same time independent property F{y). The case of ordinary statistical 
mechanics corresponds to 5 = 1/2 with F{y) being a Gaussian function of y. The case 
under study in this paper is not ordinary because F{y) is a Levy distribution and the 
scaling parameter 6 is given by 

This thermodynamic condition, however, is the asymptotic limit of a very slow transition, 
driven by the correlation function of Eq. (||), whose lifetime, in the case here under 

study, with 2 < /i < 3, is infinite. For a further discussion of this extremely slow transition 
from dynamics to thermodynamics, the reader can consult also Ref. |jl^. Although this 
process of transition to thermodynamics, or aging, is extremely slow, with virtually an 
infinite lifetime, it is not easy to detect by mean of the current techniques of analysis of 
time series, which rest on the observation of a single sequence (see for instance Ref. [0). 
Here we show that it turns out to be difficult even with the most advanced method of 
scaling detection, the method of Diffusion Entropy (DE) The first step of all these 

techniques, including the DE method, is based on deriving the random walkers required 
by Bob's experiments from a single sequence. We use Bob's algorithm to create a single, 
virtually infinite, trajectory. Then we span infinitely many portions of this sequence moving 
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along the sequence a window of width t. The selected walking trajectories are shifted in 
such a way as to make them start at a; = at t = 0. These trajectories tends to depart the 
ones from the others, as an effect of their partially random character, and this spreading 
is described by p{x,t), evaluated numerically with a proper partition of the x-axis. At 
this stage, rather than evaluating the variance, which leads to a wrong scaling in the Non- 
Gaussian case |T^, we evaluate the Shannon entropy 

/ + 00 
p(x,t)log(p(x,t)). (13) 
-oo 



In the case when the scaling condition of Eq. (|ri|) applies, by plugging Eq. (|TT]) into Eq. 



(p!3|) we get 

S{t) = A + 5\ogit), (14) 

with A being a constant whose explicit expression is of no interest here. Fig. 1 proves that 
the DE method is successful in detecting the predominant asymptotic scaling of Eq. (p!2[). 
This is a flattering result, since the methods based on the variance measurement cannot 
reveal this Levy scaling. However, it is not easy to relate Fig. 1 to the slow transition 
process described by Eq. (|TUp |T^ . 



We can go beyond the limitations of the DE method of analysis detecting the multifractal 
properties generated by the aging process itself. To reveal the emergence of these multifractal 



properties we rest on Eq. ([T0|) and on a procedure reminiscent of that of Nakao |T^. This 
author proved in fact that a truncated Levy process yield bifractal properties. It has to be 
pointed out, however, that there is a deep difference between the case considered by Nakao 
and the dynamic truncation of this letter. The effect observed by Nakao has to do with 
the ultraslow convergence to Gaussian statistics discussed some years ago by Mantegna and 
Stanley |T^, observed independently by the authors of Ref. [|T9|. The effect here discussed 
is instead an ultraslow convergence to Levy statistics. Moreover, as we have seen, our 
transition process has an infinite lifetime, while the lifetime of that of Refs. |]IB|,|TU|, although 
impressively large, is finite. 



To study the multifractal properties associated to the dynamic approach to Levy pro- 
cesses, we follow the prescriptions of Refs. ||2yj21| . We study theoretically and numerically 
the fractional moment which is expected to yield 

r+co 

(IxH = / p{x,t)\x\''dx ^t^^. (15) 



The power index C,q plays a critical role. According to the theory of Refs. pO| , ^ , C,q as a 



function of q, would be a straight line if the monofractal condition applied, while its deviation 
from a straight line signals the occurrence of multifractal properties. Using the theoretical 
prediction of Eq. ( p!0D it is straigforward to predict that 



U = 5q, (16) 

for g < /i — 1, and 

= g - + 2, (17) 

for q > fi—1. These theoretical predictions are very satisfactorily supported by the numerical 
results, as proved by Fig. 2. 

In conclusion, we define aging the process of transition from dynamics to thermodynamics 



when this transition has an infinite lifetime [22,231. If it had a finite lifetime, it would be 



possible to make an observation at times so large as to make the scaling condition of Eq. (|TT]) 
virtually exact. This condition, in turn, plugged into Eq. (|15]), would make Eq. ([T6| ) valid 
for all g's, and not not only for the small ones. This means that no aging is equivalent to the 
condition of monofractal statistics. If, on the contrary, aging exists, the diffusion process 
becomes bifractal. This is what we mean by aging induced multifractality. We propose 
aging as a paradigm for the living state of matter. This has to do with the widely accepted 



conviction (see, for instance, Ref. p3]) that life is a balance of order and randomness. This 



paper shows that aging, as we do mean it, is an attractive expression of this balance. 
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FIG. 1. Diffusion Entropy as a function of time t. Dots: numerical evaluation of S{t) for Bob's 
experiment with /x = 2.5, T !=a 1.0. Bob's walkers are derived from a single trajectory of total 
lenght 25,918,673 time units (see the text for details). The solid line represents a best fit, in the 
asymptotic regime, of the theoretical prediction S{t) = 2/31og(i) + const. 
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FIG. 2. Multifractal index as a function of q. Dots: numerical evaluation of for Bob's 
experiment (dynamical truncated Levy process) using the same data as those of Fig.l. Dashed 
line: = 5q, with 5 = 2/3 according to (|l^); dotted line: £,q = q — (3, with [3 = ^ — 2 = 0.5. 
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